Solitary Waves of Planar Ferromagnets and the Breakdown of the Spin-Polarized 

Quantum Hall Effect 
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A branch of uniformly-propagating solitary waves of planar ferromagnets is identified. The en- 
ergy dispersion and structures of the solitary waves are determined for an isotropic ferromagnet 
as functions of a conserved momentum. With increasing momentum, their structure undergoes a 
transition from a form ressembling a droplet of spin-waves to a Skyrmion/anti-Skyrmion pair. An 
instability to the formation of these solitary waves is shown to provide a mechanism for the electric 
field-induced breakdown of the spin-polarized quantum Hall effect. 
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PACS numbers: 75.10.Hk, 11.27.+d, 73.40.Hm 

Continuum descriptions of the static and dynamic 
properties of ferromagnets give rise to many forms of non- 
linear structure Q. These structures depend strongly 
on the symmetries and dimensionality of the ferromag- 
net. In two dimensions, the existence of a topolog- 
ical classification of spin-configurations (2j has impor- 
tant consequences. In particular, the topological index 
plays a key role in the construction of a sequence of 
static "Skyrmion" solitons of an isotropic ferromagnet 
||. These solitons have been of recent interest in quan- 
tum Hall systems, for which the low-energy charged ex- 
citations can be described as Skyrmions Q . 

The topological density is of particular importance for 
the dynamical properties of planar ferromagnets. It has 
been established that the dynamics are contrained not 
only by the existence of the topological index, but also 
by conservation laws based on moments of the topolog- 
ical density ||. Here, we use the properties of a con- 
served momentum, defined in terms of the dipole mo- 
ment of the topological density, to identify a branch of 
uniformly-propagating solitary waves of planar ferromag- 
nets. At large momentum the solitary waves appear as 
Skyrmion/anti-Skyrmion pairs, which move at a constant 
velocity perpendicular to the line separating their cen- 
tres. To the best of our knowledge, these are the first 
propagating solitary waves of a planar ferromagnet to be 
identified. We further show that an instability to the cre- 
ation of these solitary waves leads to a breakdown of the 
spin-polarized quantum Hall effect, akin to the destruc- 
tion of superfluid flow by vortex-ring nucleation || . 

We start from the dissipationless Landau-Lifshitz 
equation, which provides a semiclassical description of 
the dynamics of the local magnetization of a ferromagnet, 
as represented by the 3-component unit vector rij(r, t) R] 
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tion convention is assumed throughout, the labels i, j,k 
running over the three spin-components and a, j3 over the 
two spatial dimensions.) 

The total energy, E, and the modulus, nf, are con- 
served by Eq.(Q). There are several other conserved 
quantities that are important in what follows: the topo- 
logical index, Q, the linear momentum, P a , and the total 
number of spin- reversals, N, defined by ^ 



Q = J d 2 r q(r), 
P a ee 4np,Je afj J d 2 r r p q{r), 
N ee pj/h ( d 2 r (1 - n z ), 



p is the density of magnetic moments, each of angular mo- 
mentum J, and E(ni,W a nj) is an energy functional of 
the magnetization and its spatial derivatives. (Summa- 
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(3) 
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where q(r) ee -^€ijk£ a pniV ' a rijV pnk is the topological 
density. 

We will focus on configurations for which the magneti- 
zation tends to n z — 1 at spatial infinity. In this case, the 
topological index Q is an integer J^H, that is conserved 
under smooth deformations of rii(r). Q is therefore a 
constant of any regular dynamics. The conservation of 
P a and N under Eq. (Q) requires certain restrictions on 
the terms in the energy functional. It is sufficient that 
the energy functional describe a translationally-invariant 
ferromagnet, with at most uniaxial anisotropy (chosen 
along the i-axis) , and that magnetic dipolar interactions 
can be neglected These conditions are met by the 
energy functionals considered below. 

Essential to our work is the definition of linear mo- 
mentum (^) in terms of the dipole moment of the topo- 
logical density. This form avoids ambiguities inherent in 
other definitions [Q (see Ref. ||), that would prove im- 
portant for the configurations studied below. We use the 
conservation of P a , N and Q to establish the existence 
of a branch of solitary waves, following a method that 
has been used to construct solitons of the 1-dimcnsional 
Heisenberg ferromagnet M. 
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Consider the spin-configurations that minimize the en- 
ergy E for given values of P a , N, and Q. By combining 
the variational equations with Eq. ([!]) one finds that these 
spin-configurations have the time evolution 
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where v a and u> are Lagrange multipliers introduced to 
enforce the constaints on P a and N. This equation 
shows that the extremal spin-configurations are travelling 
waves, that translate at a uniform velocity v a while the 
magnetization precesses around the z axis at an angular 
frequency u>. Defining an "energy dispersion" Eq(P , N) 
by the minimum energy at given P a N and Q, the vari- 
ational equations further require 
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Whilst one can imagine searching for such travelling 
waves within each topological subspace, note that the 
momentum (U) transforms as P a — > P a — QR a un- 
der rii(r a ) — > rii(r a + R a ). Consequently, for a ferro- 
magnet that is translationally-invariant, Eq^. (P,N) — 
Eq_£ O (0, N) and all modes have vanishing velocity (^). 
For the Q = sector, however, one does expect the energy 
of a translationally-invariant ferromagnet to depend on 
the momentum. In the following, we will study the prop- 
erties of a translationally-invariant ferromagnet within 
this Q = subspace (henceforth we drop the subscript 
on Eq). Amongst the extremal spin-configurations that 
determine E* (P a , N) , we will find configurations that are 
localized in space and have non-zero v a and lo. These are 
therefore uniformly-propagating solitary waves Q of the 
ferromagnet. 

Specifically, we study an isotropic ferromagnet, for 
which the energy functional for long-wavelength distor- 
tions is the 0(3) non-linear <r-model 



V Q n,] = ip,/" d 2 r (V Q n,) 2 



(7) 



characterized by the spin-stiffness p s . Our results also 
apply to extensions of this model involving additional 
terms that are functions of Q, N and P a (these sim- 
ply shift the velocities and precession frequencies of the 
modes we construct). In the following, we will work with 
units such that p s = pj = K = 1. 

We have investigated the energy dispersion E* (P a , N) 
of the non- linear a-model numerically, using a lattice dis- 
cretization of the continuum ferromagnet (details will be 
presented elsewhere ||). The results we present are for 
a square region of 249 x 249 spins, bounded by spins 
fixed to n z = 1. The parameters are chosen such that 
the results well represent the continuum theory (^) de- 
fined within a square region (of side L, say). The scale- 
invariance of such a theory requires that E* (P, N, L) = 
E*(PX, NX 2 ,LX) for all scale-factors A, so the energy can 



depend only on the "scaled momentum" p = Pj y/N, and 
the "boundary parameter" r] = N/L 2 . The properties of 
an infinite system are determined by considering the limit 
77 -> 0. 
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FIG. 1. The energy dispersion of the isotropic ferromagnet, 
for n = 0.032 (solid line) and n = 0.064 (crosses). The en- 
ergies, for an infinite system, of free spin-waves (dashed line) 
and a non-interacting Skyrmion/anti- Skyrmion pair (dotted 
line) are shown for comparison. Inset: the same quantities as 
functions of 1/p , indicating the three regimes of the n — > 
limit. 



Figure [l] shows results for the minimum energy as a 
function of p at two values of the boundary parameter 
n. Although weak dependences on rj do remain, it is 
clear that the energy dispersion E* (p, 77) is tending to a 
smooth function as 77 — > 0. The limiting configurations 
will describe travelling waves with non-zero velocity and 
precession frequency (p). 
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FIG. 2. Variation of the square of the polar angle at the 
centre of the system with scaled momentum. An extrapola- 
tion, using the filled circles, suggests that as n — > 0, 9 C would 
vanish for p < 3.4. Also shown are the results for an ansatz 
described in the text. 
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A careful study of the spin-configurations shows that 
for 77—5-0, there is a transition from spatially-extended 
configurations at small p, to spatially-localized configu- 
rations at large p. This is illustrated in Figure which 
shows the (square of the) polar angle at the centre of 
the system 9 C as a function of p for several values of 
77. For p > 3.5, 9 C is weakly-dependent on 77, and can 
be assumed to represent well an unbounded system. For 
p < 3.5, 9 2 c decreases steadily down to the smallest values 
of 77 studied, showing that the configurations remain sen- 
sitive to the presence of the boundaries. An extrapolation 
indicates that 9 2 would vanish at a scaled momentum of 
p* = 3.4 (the dashed line in Fig. |). 

Within the assumption that 9 C does vanish as rj — > 
for p < p* , the properties in this regime can be un- 
derstood within a linearized continuum theory in which 
all quantities are expanded to second-order in the po- 
lar angle 9{r). The linear theory has extremal solutions 
9{r) = 9 cos(ttx/L) cos(ny/L), </>(r) = ky, which lead to 
an energy dispersion EswiPiV) = P 2 + 2ir 2 r]. The ener- 
gies of these configurations coincide with the energy of N 
non-interacting spin- waves confined to a square of side L. 
We refer to them as "free spin- wave" states. The numer- 
ical results for p < p* are consistent with the minimal- 
energy spin-configurations of the infinite system being 
free spin-wave states: the spin-configurations ressemble 
those of the linearized continuum theory and their en- 
ergies approach the free spin- wave energies of an infinite 
system continuously as 77 — > (see Fig. ^). 

In contrast, the spin-configurations above the transi- 
tion, p > p* , cannot be described within the linearized 
theory. In this regime, 9 C tends to a finite value as 77 — > 0, 
and the spin-configurations are spatially-localized. These 
configurations constitute the branch of propagating soli- 
tary waves that forms the main subject of the present 
paper. 
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FIG. 3. Spin-configurations of the solitary waves at (a) 
p — 3.64, and (b) p — 12.4. The vectors show the projection 
of the spin onto the a;y-plane for those spins with n z < 0.98 
(for clarity only one spin in six is plotted). The direction of 
the momentum and velocity of the solitary waves is indicated. 

Typical spin-configurations of these solitary waves are 
shown in Fig. ||. For small momenta (yet larger than 
p* ) the configurations ressemble the free-spin- wave states 
described above, but are localized in space [Fig. ||(a)]. 



We refer to these configurations as "spin- wave droplets" . 
For large momenta, two vortices of opposite sign ap- 
pear in the projection of the spins on the xy plane 
[Fig. pKb)]. In fact, such spin-configurations contain no 
singularity, since the vortices occur around points at 
which n z = —1, as occurs in (anti-)Skyrmion configu- 
rations ||. Accordingly, we refer to these configurations 
as "Skyrmion/anti-Skyrmion pairs". The transition be- 
tween spin- wave droplets and Skyrmion/anti-Skyrmion 
pairs is found to occur at p = 4.8 (see Fig. ||). Note 
that this transition is effected by continuous deforma- 
tions of the spin-configuration. The spin-wave droplet 
and Skyrmion/anti-Skyrmion pair configurations there- 
fore form part of the same branch of solitary waves; the 
differing names we have assigned indicate only the differ- 
ences in the forms of the spin-configuration. 

As the scaled momentum increases beyond p = 4.8, 
the separation between the Skyrmion and anti-Skyrmion 
grows continuously and the energy approaches the en- 
ergy, 87r, of a non-interacting Skyrmion/anti-Skyrmion 
pair H . Asymptotic forms that have been proposed for a 
widely-separated Skyrmion/anti-Skyrmion pair have 
a number of spin-reversals N that diverges logarithmi- 
cally with system size, and so correspond to the scaled 
momentum p = 0. However, deformed versions of these 
configurations can be expected to lead to asymptotic in- 
teractions of the form ~ l/[p 2 \og(p)] ||. We find an 
approximately 1/p 2 dependence (see the inset to Fig. [I]). 
Logarithmic corrections cannot be studied, owing to 
finite-size effects that appear once the Skyrmion/anti- 
Skyrmion separation ~ P/Ait becomes comparable to the 
system size L (i.e. when p > 4 ir/^/rj). 

Previous numerical studies UM of the time- 
development of a spin-configuration ressembling Fig. ||(b) 
under a Landau-Lifshitz dynamics showed evidence for 
a uniform drift parallel to the momentum. The present 
work establishes that solitary waves of this form exist for 
the non-linear cr-model, and determines their structures 
and velocities. 

Having described the properties of the free spin-waves 
and solitary waves, we return briefly to discuss the tran- 
sition at p* = 3.4 that was inferred by the extrapolation 
shown in Fig. ||. Further evidence that a transition oc- 
curs is obtained by making a radially-symmetric ansatz 
for the spin-wave droplet states: 9(r) = /(H), <j> = ky. 
The non- linear equation for /(r) that arises from ex- 
tremizing the energy has been studied in another context 
p3|. In the present context, the results indicate that lo- 
calized solutions are lost below p* ad = 3.4 (see Fig. ||). 
The existence of this transition lends weight to the idea 
that a transition occurs in the unrestricted case. Fur- 
thermore, the coincidence of the two transition points, 
p* = p* ad , suggests that the configuration at the transi- 
tion for the unrestricted problem is of the form assumed 
in the radially-symmetric ansatz. 

We now turn to discuss the application of these re- 
sults to quantum Hall ferromagnets (QHFs). These are 
the incompressible states responsible for the quantized 
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Hall effect of two-dimensional electron systems at certain 
filling fractions [fl3|]4"|. The long- wavelength dynamics 
of QHFs are described by the Landau-Lifshitz equation 
[§,|h|, with the Zeeman energy, (gpBB)N , and the non- 
linear CT-model term (f^) as the leadin g te rms in a gradient 
expansion of the energy functional |lq| . We will adopt 
this level of approximation, allowing the results presented 
above to be applied directly. Higher-order terms in the 
gradient expansion (e.g. long-range Coulomb forces Q|) 
will not affect the qualitative discussion [||. 

Our results provide a semiclassical description of the 
neutral excitations of a QHF involving N spin-reversals 
with total momentum P (this description is accurate in 
the limit N ^> I). The nature of the low-energy ex- 
citation depends on the value of the scaled momentum 
p: for p < p* it is a collection of N spin-waves, each of 
momentum P/N; for p > p* it is a solitary wave, ressem- 
bling a spin- wave droplet or Skyrmion/anti-Skyrmion 
pair. These solitary waves are highly reminiscent of the 
"magnetic exciton" picture of a single spin- wave |l6| : two 
oppositely-charged particles move at a constant velocity 
perpendicular to the line joining their centres. Indeed, 
just as Skyrmions are the large-spin versions of the quasi- 
particles of a QHF, these dynamical solitary waves are 
the large-spin analogues of the spin-waves. 

Quantum effects, not accessible within the semiclas- 
sical description, could lead to certain differences. If 
the quasiparticle gap is less than the Skyrmion gap for 
vanishing Zeeman energy JlTj, the lowest-energy states 
at large scaled momentum will not be Skyrmion/anti- 
Skyrmion pairs; they will be widely-separated quasielec- 
tron/quasihole pairs. Also, if two spin-waves are weakly- 
bound at any non-zero total momentum then con- 
figurations at all values of p will be spin- wave droplets. 
The transition found at p = p* will then be a crossover, 
becoming a true transition only in the limit N — > oo. 

Experimentally, one way in which the solitary waves 
can be accessed is through the application of an electric 
field £ a to the QHF. For weak electric fields a quan- 
tized Hall current will flow with exponentially-small dis- 
sipation (we assume the limit of vanishing temperature). 
However, there exists a threshold value of the electric 
field |l9|, £t = ji\/167rp s gfiBB, above which the sys- 
tem can reduce its energy by creating single spin-waves 
of small momentum. The spontaneous creation of spin- 
waves is expected to initiate the breakdown of the quan- 
tum Hall effect Our results indicate that, as the 
number of spin- waves increases, one can expect spin-wave 
droplets and Skyrmion/anti-Skyrmion pairs to form. The 
appearance of Skyrmions was anticipated in a footnote 
to Ref. jnj . Our results provide a detailed description of 
the non-linear structures that can develop. 

Firstly, we can determine the number of spin-waves 
involved in the formation of a spin- wave droplet. For a 
QHF placed in a field equal to the threshold field £t (such 
that single spin- waves of momentum y/ gpsB/4:Trp s fi/£ 
are formed), it is energetically favourable for spin- waves 



to collapse to a spin-wave droplet when their number 
reaches N c — p* 2 p s / (gpBB). 

Secondly, our results show that the spin-wave droplet 
that forms at this point is unstable to increasing its mo- 
mentum (since the solitary-wave group velocity is maxi- 
mal at p — p* ) . If energy can be lost from the electron gas 
to the lattice, the droplet will be torn apart in the electric 
field, forming a Skyrmion/anti-Skyrmion pair of increas- 
ing separation (these charges will eventually reach the 
voltage contacts). By the continual nucleation of spin- 
wave droplets, and their relaxation to widely-separated 
Skyrmion/anti-Skyrmion pairs, a d.c. current can be car- 
ried along the electric field. The quantized Hall effect is 
thereby destroyed, much as superfluid flow is destroyed 
by the nucleation and growth of vortex rings || (note, 
however, that for £ > £t the spin-wave droplet can nu- 
cleate without thermal fluctuations). 

This breakdown mechanism could most effectively be 
investigated in GaAs samples under hydrostatic pressure, 
for which the g-factor, and thus £t, can be made to 
vanish |^(J. In realistic systems, disorder will affect the 
breakdown mechanism strongly. Most importantly, the 
localization of the Skyrmion and anti-Skyrmion will pre- 
vent the run-away separation expected for a clean sys- 
tem, and charge transport will proceed by the conden- 
sation of spin-wave droplets in the presence of pinned 
Skyrmion/anti-Skyrmion pairs. 

I am grateful to Fabio Pistolesi, Tim Ziman and John 
Chalker for helpful discussions and advice. 
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